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Abstract 

We calculate the cross-section for the thermal n + p ^ d + j process in chi- 
ral perturbation theory to next-to-next-to-leading order using heavy-fermion 
formalism. The exchange current correction is found to be (4.5 it 0.3) % in 
amplitude and the chiral perturbation at one-loop order gives the cross sec- 
tion a^^ = (334 lb 2) mb which is in agreement with the experimental value 
(334.2 lb 0.5) mb. Together with the axial charge transitions, this provides a 
strong support for the power of chiral Lagrangians for nuclear physics. 
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One of the corner-stones of nuclear physics is the successful explanation in terms of 
exchange currents given two decades ago by Riska and Brown [Q of the ~ 10% discrepancy 
between the experimental cross-section and the theoretical impulse approximation prediction 
for the process 

n + p^ d + 'y (1) 

at threshold. Riska and Brown computed, using a realistic hard-core wave function for 
the deuteron, the two one-pion-exchange diagrams initially suggested in 1947 by Villars 
plus the uj and A resonance diagrams. That the dominant contributions to electroweak 
exchange currents could be gotten from current-algebra low-energy theorems was suggested 
by Chemtob and Rho who gave a systematic rule for organizing the leading exchange- 
current diagrams effective at low energy and momentum. Although suspected since the 
Yukawa force was introduced, the work of Riska and Brown was the first unequivocal evidence 
for the role of mesons, in particular that of pions, in nuclear interactions. In this Letter, 
we show that the terms considered by Riska and Brown are a (main) part of the terms 
that figure in chiral perturbation theory to next-to-next-to-leading (N^L) order and that 
when completed by the rest of the N^L order terms, chiral perturbation theory scores an 
impressive success in nuclei. 

In the modern understanding of QCD, it is the spontaneous breaking of chiral symmetry 
associated with the light quarks that predominantly governs the structure of low-energy 
hadrons as well as the forces mediating between them. In fact, the full content of the 
gauge theory of strong interactions, QCD, can be expressed at low energy by a systematic 
chiral expansion starting with effective chiral Lagrangians @]. Stated more strongly, such 
an approach, known as chiral perturbation theory (xPT), while reproducing the current 
algebra, is now considered to be exactly equivalent to QCD in long wavelength regime . 
Our paper reports the first quantitative chiral perturbation calculation of the fundamental 
nuclear process and shows that chiral symmetry is indeed a powerful guiding principle 
in nuclear dynamics, confirming the work of Riska and Brown ||l[ and the conjecture of 
Kubodera, Delorme and Rho 0. 

Two recent developments provide a strong motivation for this work. The first is the work 
of Weinberg and Ordofiez, Ray and van Kolck on understanding nuclear forces from 
chiral Lagrangians. The second is the explanation by the present authors [^] of the enhanced 
axial-charge transitions in heavy nuclei in terms of exchange currents in chiral perturbation 
theory treated to the same chiral order as for nuclear forces. 

In both cases cited above, one is limited to long wavelength processes, with the typical 
energy /momentum scale Q much less than the chiral symmetry scale ~ my ~ 1 GeV. 
This is because xPT is an expansion in Q/A.^ and its practical value lies where Q/A^ ^ 1- 
This entails a subtlety in applying xPT to nuclear processes as we shall now specify, a 
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feature absent in such "elementary processes" as vrvr or vrA^ scattering for which much work 
with impressive success has been done [[lO| , p!T| . 

The precise way xPT can be apphed in nuclear physics was explained in |^. Here we 
sharpen the key arguments to bring home our thesis. For nuclear physics, where baryons 
as well as mesons are involved, the chiral expansion is made in heavy-fermion formalism 
]T^,[7|] which allows a systematic expansion in derivative on pion fields, d^^/A^ as well as on 
baryons fields, Bb/X with X = A^^ or (baryon mass) and in ttIt^/K^ where m^r is the 
pseudo-Goldstone boson mass. The expansion is organized by the power v in Q*^, given by 
Weinberg |^, 

V = A-Nn-2C + 2L + Y,A, (2) 

i 

where Aj = (ij + |nj — 2 and Nn is the number of nucleons involved, L the number of 
loops, di [rii) the number of derivatives or powers of pion mass (nucleon lines) that enter 
into the ith vertex and C the number of separated pieces of the Feynman graphs. In the 
presence of external fields {i.e, electroweak currents) which is what we want to study here, 
chiral invariance requires [jl3[ that Aj > —1. The xPT in nuclear systems amounts then 
to compute Feynman diagrams involving external fields in the increasing power u embedded 
inside the most general process describing the transition from the initial nuclear state to 
the final nuclear state with interactions taking place before and after the current insertion. 
This is essentially what was first suggested in but on a somewhat ad hoc basis. What 
this means is that we are to take the most realistic nuclear wave functions and calculate the 
transition matrix elements with the xPT graphs computed to the maximum possible order 
of chiral expansion. It is in this sense that the "counter terms" in the chiral expansion can 
be fixed from experiments. This point has been emphasized also by Weinberg in his 
discussion of nuclear forces, in particular many-body forces. An important consequence of 
this strategy is then that only the current operators obtainable by xPT are to be kept. This 
implies that (1) short-wavelength effects are to be "filtered out" and (2) n-body currents with 
n> 2 are suppressed in the same sense that ra-body forces are suppressed 0. Of particular 
importance of the first implication is that when the currents are put in coordinate space, 
shorter-range interactions for < where is the hard core radius cannot contribute 
in xPT . If the strategy is correct, the result should not sensitively depend on the precise 
value of Tc, within the relevant range for application in nuclei to the order considered, say, 
A;^^ < rc < (2m^)-\ Note that this is roughly the range that has been successfully described 
in xPT for the NN potential [g. 

We now focus on the process @) at thermal energy. The relevant operator is the isovector 
magnetic moment operator which we shall denote /i. The vector current for the two-nucleon 
process consists of the one-body current J^x (called "impulse current" ) and the two-body 
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current Jjt^-^ (called "exchange current"). Although the time part of the leading order one- 
body current has z/ = — 3 according to the counting rule (j^), the space part - which is 
relevant to the /i - has u = —2, suppressed by a factor of order 0{Q/m]\f) for exactly the 
same reason as for the suppression of the time component of the single-particle axial-current 
0]. The leading order two-body current has u = —1 with one- loop corrections entering at 
u = +1. Thus from the point of view of chiral expansion, the two-body current at one- loop 
order corresponds to next-to-next-to-leading (N^L) order for the magnetic moment operator. 
This is the order that was computed in the case of the nuclear axial-charge transitions studied 
in Ref. and that will be adopted for the vector current matrix element for the process 

(B). 

To the N^L order, the one-body current is the usual impulse isovector magnetic moment 
operator with the renormalized vertex i.e., = f^v J2i '^i where is the proton 
mass and fiy = fip — fi^n — 4.70589. Now to the next-to-leading order {i.e., 0{Q) relative 
to the impulse), the two two-body currents Fig.l(al) and (a2) contribute. With the renor- 
malizations of the vertices and of the nucleon lines, they give the "tree two-body current" 
involving the renormalized (measured) constants qa = 1.257 and /,r = 93 MeV. The result- 
ing operator will be denoted by fi(2)tree- We should mention that this is the same as what 
one gets from the corresponding "pair" and "pionic" currents given in with the renormal- 
ized ttA^A^ coupling constant g^NN replaced by rnjsgA/ f-K through the Goldberger-Treiman 
relation. 

Going to the N^L order {i.e., 0{Q^) relative to the impulse), we encounter two classes of 
two-body currents: one-pion exchange with one-loop radiative corrections at the vertices and 
two-pion exchange. Part of the former class of diagrams renormalize the soft-pion exchange 
current to give the tree two-body current mentioned above. Besides this trivial correction, 
there is a nontrivial contribution coming from finite counter terms in the irVNN vertex in 
Fig.l(al), 
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+ 



' B (v'^v'' - g^^ - \s^S'') {5ab - ^ j r^„5" 5 + h.c. 



4/2 n - 3 



(3) 



where T = f = d^T, - d,T ^ + [F^, T,] and = f F^ (A^ = f A^) is a vector 
(axial-vector) covariant combination of pion fields which contains either one derivative or 
an external gauge field. It can be shown that the counter-term constants c and c' can be 
obtained by saturating with the resonances uj and A respectively, 

g^mi 

In 8n^gA{mi - mi) 

CA - ^ C = ^ 0.1667 (4) 

9/^ 9gA{mA - mNjiriN 
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where gu, is determined from the — 717 decay, guj = 17.55, and the A^A transition magnetic 
moment /it and the vrA^A couphng C come from the fit to the A properties as explained in 
TM, C = —1.73 and /iy = —7.7 ± 0.5. The resulting one-pion exchange two-body current 



is given by the two graphs Fig.l(a3) and (a4). We shall denote the corresponding magnetic 
moment operator by ^{2)Itv While these are again identical to what was obtained in [Q], their 
interpretation is entirely novel and significant: chiral symmetry tells us that they constitute 
the complete (isovector) one-pion exchange current corrections to /i(2)tree- In particular, there 
are no finite chiral loop corrections to the one-pion- exchange operator. This is in contrast to 
the case of axial-charge transitions 0. It seems reasonable therefore to lump all four terms 
of Fig. la together and call them "generalized tree operators." 

There are numerous diagrams of genuine loop character that can contribute in general 
kinematics to the two-pion-exchange two-body current. Fortunately things simplify dras- 
tically for the process (|l]) in heavy-fermion formalism, with only four graphs (bl), (b2), 
(b3) and (b4) of Fig.l non- vanishing: Two-body currents involving four- Fermi interactions 
are zero-ranged in coordinate space and together with all other zero-range terms give zero 
contributions due to the short-range cut-off (or correlation) as described in detail in 0. 
Now these four graphs - which appear at the same chiral order as the above "generalized 
tree operators" - have hitherto been unaccounted for in this form in the previous studies of 
exchange currents. (It is possible however that part or all of this may have been included 
in heavy- meson exchange graphs in phenomenological approaches.) After divergences are 
removed (we use dimensional regularization) , the resulting magnetic moment operator takes 
the form 

/^(2)2. = (^T^-) - r(-)) j^dttVW^y,{2tx^) 

^^l9A nj2^ix) ( y,i^2tX.) -}-y,i^2tX. 





where x^ = m^r, yo{x) = yi{x) = (1 + x)yoix) and y2ix) = (l + ^ + yoix) . 

The spin-isospin operators are defined by T^®"* = (ri T2)^(cri (T2) and Tj,®-* = (ri 



r2Y 



rf ■ ((Ti (J2) — |(cri (T2) for = — and x. 



The total magnetic moment operator to the N^L order then is 

+ /^(2)l7r + /^(2)27r- (6) 

Let the corresponding matrix elements for the capture process be denoted by Mj. The we 
are specifically interested in the ratios 6tree = ^{2)tree/^(i), ^itt = ^(2)i7r/^(i)5 S{2)2tt = 
M(2)2tt / M^i) and 823 which is the sum. 
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We now describe the numerical results. For this, we should pick the most realistic two- 
nucleon wave functions. We take the Argonne potential vis recently constructed by Wiringa, 
Stoks and Schiavilla This potential is fit to 1787 pp and 2514 np scattering data in 

the range 0-350 MeV with an excellent of 1-09 and gives the deuteron properties - the 
asymptotic S-state normalization, As, the D/S ratio, r], and the deuteron radius, da - 
close to the experimental values. Electromagnetic properties also come out well, modulo 
exchange- current and relativistic corrections. 

We use the physical values for masses and constants that appear in the theory. There 
are no unknown parameters except one that has to do with nuclear interactions at short 
distance, namely short-range cutoff Vc- Figure 2 summarizes the results of the calculation 
for a wide range of Vc, < Tc < 0.7 fm compared with the experiment |l^. We describe 
very briefly how this result comes about. Details will be given in a later publication . 

The capture cross section is proportional to where a<j is the ^5*0 np scattering length. 
The Argonne vis predicts al^ = —23.732 fm in excellent agreement with the experimental 
value af^ = —23.739 ± 0.008 fm. The single-particle matrix element with this potential 
gives the impulse approximation cross section (given by of Eq.(Q)) aimp = 305.6 mb, 
about 9.6% less than the experimental value aexp = 334.2 ±0.5 mb. It has no Vc dependence. 
The two-body matrix elements computed with the same wave function does depend on 
but very weakly. The resulting total cross section is plotted in Fig. 2a for the relevant range 
of Tc- (The Tc dependence is incorporated in the radial integral by multiplying the integrand 
by ^^(r — Tc).) Figure 2b shows the contribution of each term in terms of the ratios 5(2). We 
note that the "generalized tree" contributions dominate to the N^L order, with only a small 
correction (less than 0.6 % of the single-particle matrix element) coming from the genuine 
one-loop correction. This agrees with the "chiral filter" mechanism seen in the axial-charge 
transitions and confirms the conjecture made in The intrinsic uncertainty associated 
with short-distance physics notwithstanding, the theoretical prediction (t^pt = 334 ± 3 mb 
(where the theoretical error bar represents the dependence on the hard-core cut-off) is in 
remarkable agreement with the experiment, say, within less than 1% ! 

In conclusion, we discuss the meaning of the short- distance cut-off in xPT . As dis- 
cussed in 0, the loop terms contain zero-range operators in coordinate space. In addition, 
four-Fermi counter terms in the chiral Lagrangian with unknown constants are also zero- 
ranged. At higher chiral order, increasingly shorter-ranged operators would enter together 
with the zero-ranged ones. Now if we were able to compute nuclear interactions to all or- 
ders in chiral perturbation theory, the delta functions in the current would be naturally 
regularized and would cause no problem. Such a calculation of course is an impossible feat. 
The practical application of chiral perturbation theory is limited to low orders in the chiral 
expansion, so a cut-off would be needed to screen the interactions shorter-ranged than acces- 
sible by the chiral expansion adopted. Clearly such a calculation would be meaningful only 
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if the dependence on the cut-off were weak. Our calculation here and also the one in |^ meet 
that criterion. In the present work, we find Vc — 0.5 fm at which the theoretical prediction 
agrees exactly with the experiment indicating that at one-loop order, chiral perturbation 
theory is meaningful for processes taking place at internucleon distances r > 0.5fm. We 
conjecture that to the next order, that is 0{Q^) relative to the leading tree graphs Fig.l(al) 
and (a2), the appropriate cut-off would be Vc ~ l/w^t^ ~ 0.3 fm. It would be a challenge to 
quantify this heuristic reasoning. 
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FIGURES 
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FIG. 1. The Feynman graphs contributing to the two-body vector current for the process (\\^: 

(a) the "generahzed tree graphs"; (6) two-pion exchange graphs. The current is depicted by the 

wiggly line, the pion by the broken line and the nucleon by the sohd hne. One-loop graphs figuring 

in ttVNN vertex are entirely saturated at low photon energy by the two resonance-exchange-tree 

graphs ((a3) and (a4)). 
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FIG. 2. Total capture cross section cjcap (top) and 5's (bottom) vs. the cut-off rc- The solid 
line represents the total contributions and the experimental values are given by the shaded band 
indicating the error bar. The dotted line gives (Jtreej the dashed line 6tree + ^i^, the dot-dashed line 
(5tree + ^in = ^trec + ^^i^ + ^itt and the solid line the total ratio, 52b- 
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